clear;
clc;
a=sym('a%d',[1,4]);
H=hankel(a);


H1=[H eye(4)];
H2=rref(H1),
H3=H2(:,5:8)